Dynamics of entanglement in two coupled qubits 
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Considerable progress has been made in understanding how dissipation drives the state of a single 
qubit from a quantum-mechanical superposition to a classical mixed state. This paper uses a Bloch- 
Redfield approach to study how, in a system of two qubits, dissipation drives the bipartite state 
of the two qubits from an entangled state to a product state. The measure used for mixed-state 
entanglement is Wootters's formula for the entanglement of formation. For qubits that start in an 
entangled state and are then decoupled, entanglement is found to decay faster in general than the 
decoherence of single qubit states by a factor 2/ log 2 that can be obtained analytically in a limiting 
case. We show that the dynamics of entanglement for realistic parameters is different from that of 
fidelity: for some experimental parameters a bipartite state loses fidelity much more rapidly than 
it loses entanglement. These results are of some practical interest since entanglement is thought to 
be the essential "fuel" for fast quantum algorithms. 

PACS numbers: 85.25.-j, 85.25.Cp, 85.25.Dq 



I. INTRODUCTION 

Over the past few years considerable progress has 
been made in fabricating and manipulating solid state 
qubits 1 ^^*^^^. A particular implementation of in- 
terest involves superconducting elements where either 
the charge or flux states are exploited for the pur- 
poses of quantum computationifi*ii*iSiii. Other sys- 
tems that are also being actively studied include trapped 
ionsJ2ii5ii£iAIii£ ) an£ j nuc i ear moments of atoms in an 
NMR setupi2i2fli2i. One of the requirements for the pur- 
poses of quantum computation is the ability to isolate 
qubits from external noise. Unfortunately perfect isola- 
tion is never possible in practice, as the solid state qubits 
are always coupled to their environment, but much expe- 
rience has been gained in observing quantum phenomena 
in macroscopic solid state systems^. The processes of 
preparing, manipulating and measuring the state of the 
qubits also provide mechanisms for decoherence and dis- 
sipation. 

The standard theoretical approach to decoherence in a 
single qubit models the environment as a bath of oscilla- 
tors coupled to the qubi1*22: the dissipative environment 
causes a coherent superposition in the qubit to evolve to- 
ward a mixed state. Current interest in the field of quan- 
tum computation is focused on fabricating coupled qubit 
systems so as to be able to perform nontrivial computa- 
tions based on quantum algorithms. Of great importance 
in multi-qubit systems is the notion of entanglement, a 
characteristic trait of quantum mechanics, which is be- 
lieved to be the source of the exponential speed-up of 
quantum algorithms^. In this paper we study the time 
evolution of entangled states in the presence of a dissipa- 
tive environment: a fundamental question about coupled 
qubit systems is how dissipation causes entangled states 
to evolve into product states. 

For single qubits, the quantum nature of the system 
is investigated by looking at either the Rabi oscillations 



or Ramsey fringe measurement, both of which rely on 
the quantum coherence of the state in question. The 
presence of an environment leads to loss of quantum co- 
herence, or dephasing, which limits the operation of the 
qubit as a device for computation. On the other hand, for 
coupled qubits, the key quantum signature is entangle- 
ment and our primary motivation is to investigate the two 
qubit system with entanglement as the diagnostic. We 
find interesting temporal behavior of the entanglement of 
formation 2 - in coupled systems, and find that once cou- 
pling is turned off, entanglement decays at a rate that 
is roughly three times faster than the single qubit deco- 
herence rate. This may prove relevant for implementing 
logic gates in quantum computers. 

Entanglement is a measure of the amount of informa- 
tion of one qubit that can be obtained by making a mea- 
surement on the second of a pair of coupled qubits. To 
illustrate consider the singlet state of a two spin system: 
= (l/\/2)(|U) - lit))- Measurement of the state of 
one, completely determines the outcome of a measure- 
ment performed on the other. Thus, the notion of entan- 
glement. Such a determination would not be possible if 
the state of the two spin system could be represented as 
a product state of the two individual spins, i.e. \ipx) \ip2}' 
any product state has zero entanglement. In particular 
the (maximally entangled) singlet state cannot be repre- 
sented as a product of individual states. The two states 
of the spins can be identified with the and 1 states of 
a qubit. 

For a bipartite system in a pure quantum state, the 
von Neumann entropy of either of the pair can be used 
as a measure of entanglement, E(ip) = — Tr(p Log2 p), 
where p is the partial trace of (ip\ over either of the 
subsystems. For a mixed state the quantification is more 
complicated, because the definition of entanglement must 
distinguish between quantum entanglement and classical 
correlation. While the singlet state can still be used as 
the basic unit of entanglement, the number of singlets 
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needed to create a mixed state is different from the num- 
ber one can extract from the state. (By extracting we 
mean the number of singlets produced by local operations 
and classical communication from an entangled state ip.) 
For the purposes of quantitative calculation of entangle- 
ment dynamics, we will follow the prescription given by 
Wootters for the entanglement of formation^, reviewed 
in Section III below. 

Gate performances and coherent dynamics in coupled 
qubit systems have previously been studie d 25 ' 26 , where 
the coupling between the two is of the form cr z ® af^ 
(Ising coupling) and Oy ® af^ respectively, er- being the 
Pauli matrices describing the two-state quantum system 
i. (In this notation it is a z that couples to the dissipative 
bath.) The models we introduce in Section II are gener- 
ally similar to these, although we find that the dynamics 
to be more interesting in the case where the coupling 
has SU{2) symmetry, i.e., is of the form crW eg) a^ 2 \ as 
is likely appropriate for spin qubits. However, the pri- 
mary emphasis of our work is the study of entanglement 
rather than fidelity of quantum evolution processes of two 
qubits. We find that some operations destroy fidelity but 
preserve entanglement, while others destroy both fidelity 
and entanglement; hence the two are sensitive to different 
physical processes. 

Preservation of entanglement is essential for applica- 
tions because lost entanglement cannot be restored by 
local operations and classical communcation. Sections 
IV and V study how entanglement evolves in the models 
of Section II, for a variety of initial states and couplings. 
In the limit of zero coupling the entanglement is seen to 
decay due to the single qubit decoherence, while at fi- 
nite coupling the system can show nonmonotonic behav- 
ior where the initial entanglement is lost at intermediate 
times, but is regained as the system evolves further. This 
behavior is indicative of a transition from one entangled 
state to another caused by dissipation, which can be in- 
ferred by the time evolution of the fidelity. 



II. THE MODEL 

Each individual qubit is a two state system that we de- 
scribe in the pscudospin representation along the z axis: 
the two states are | |), \ |). For flux qubits, the bias e 
along this direction is proportional to the net flux, <E>, 
through the superconducting loop making up the qubit, 
e = 2I p (Q — ^$o)j where I p is the persistent current in 
the loop and <£>o is the flux quantum. In addition there 
exists a tunnel splitting, A, between the two states with 
currents flowing in opposite directions at the degeneracy 
point, <!> = $o/2. This is captured in the Hamiltonian 
with a term proportional to the rr-component of the pseu- 
dospin. Similar representations for the Hamiltonians for 
the NMR and ion trap realizations for qubits can be de- 
veloped where the bias is the z component of the applied 
static magnetic field and the tunneling A is the time de- 
pendent transverse magnetic field, which has components 



both in the x and y directions. 

The environment is modeled by a bath of oscillators 
which is coupled to the individual qubits. The form of 
the coupling depends on the system in question. For flux 
qubits, the noise in the flux threading the superconduct- 
ing loop is the primary source of dissipation. Hence the 
Hamiltonain includes a term that couples the oscillator 
degrees of freedom to that of the ^-component of the 
pseudospin. Similar models would apply for the NMR 
and ion trap systems if noise in the static applied field 
is the primary source of decoherence. The final piece of 
the model is the coupling between the two qubits. For 
NMR setups, there are two forms of interaction; the spin 
dipolar coupling and the through-bond interaction (or 
J-coupling). The former averages away leaving only a 
Heisenberg-type interaction between the two spins. For 
flux qubits inductive coupling leads to term of the type 
dy eg) while for other superconducting realizations 

a coupling of the form a z Cg> a z is appropriate. Here we 
study the Heisenberg and Ising cases as representative 
examples. Putting all the pieces together, we arrive at 
the following Hamiltonian (with Heisenberg coupling) 



H — H q + Hb at h + H qq + H q ^bath (1) 

Hq = -i ei a« - \^ - - U 2 a^ 



n 11 2 



Hq-bath = -^C lf f 2 (1) X 4 + i^C l( f z (2) X, 



where ej's are the biases, Aj's are the tunnel split- 
tings, J is the coupling strength and m and to char- 
acterize the bath of harmonic oscillators coupled to 
the two qubits. Individual qubits are biased such 
that ej w 5A, ; . We will further simplify to the case 
where the two qubits are identical. The eigenvalues are 
(| J - iVe 2 + 2n 2 ,-§J, | J, |J + iVe 2 + 2n 2 ), where 
e = ei + 62 and n = (Ai + A2)/V2. The correspond- 
ing eigenvectors in the singlet triplet basis (|"|"f) , (|T4) + 
IT))/V2,|U),(|Tl>-|T>)A/2) are, 
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where tan# = e/\^2n. Note that in the limit A = 0, the 
eigenvectors are the singlet and the triplets themselves. 
For purposes of our later discussion, we will continue to 
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refer to them as the singlet and triplet, as the operation 
point for individual qubits is such that A/e <C 1. We will 
also present results for the case where the interaction is 



of the form H qq = Ka 



(i) 



J2) 



III. ENTANGLEMENT 

To study entanglement, we will follow the prescription 
of Wootters^i Given a mixed state p of two quantum 
systems A and B, consider all possible ways of expressing 
the density matrix, p, as an ensemble of pure states. For 
states \tpi) and probabilities Pi, 



(3) 



The entanglement of formation is defined as the mini- 
mum, over all such ensembles, of the average entangle- 
ment of the pure states making up the ensemble 



(4) 



To avoid the process of extremization, Wootters derived a 
formula for entanglement 2 ^. Consider the density matrix 
p and its complex conjugate p* taken in the standard 
basis for a pair of spin-^ systems (||T) i Iti) > ITT) j ITT))- 
Define the spin flipped state as 

p = (cr v ® (T y )p*(a y <£> a y ). (5) 
The entanglement of formation is given by 



E{C{p)) 



i + V^c 2 



where 



h{x) — — x log 2 x ~ (1 — x) log 2 (l — x) 



C(p) = max[0, Ai - A 2 - A3 - A 4 ] 



(6) 



(7) 



and the AjS are the eigenvalues, in decreasing order, of 
the Hermitian matrix 



R = v VppVp- 



(8) 



To get a better understanding of this definition of en- 
tanglement, consider first pure states. We will work in 
the following basis (Bell basis), 



\ei) = ^(|TT) + |U)) 

lea) = ^ (ITT)- ITT)) 

|e 3 ) = i^(|U> + UT» 

|e 4 ) = ^(lU)-UT))- 



(9) 



Consider a pure state = J2i=i a i \ e i)- The von Neu- 
mann entropy of the reduced density matrix of either one 
of the qubits is given by JHl where the concurrence, C, 



is C 



E 



„2 I 28 



The same result is obtained by using 



the definition of C, via the eigenvalues of the matrix R. 
Furthermore, the trace of the matrix R is a measure of 
"distance" between the states p and p, which in turn in- 
dicates how closely the state p approximates a mixture 
of bell states22i2S. Thus the measure of entanglement 
reviewed here is a generalization of the von Neumann 
entropy for the pure case. 



IV. TIME EVOLUTION OF THE DENSITY 
MATRIX 

In the presence of dissipation the time evolution of the 
density matrix can be computed using the Bloch-Redfield 
formalism 31 . The equation of motion for the p is 3 ^ (Here 
we follow the notation in 2 ^) 



dpnm(t) 

dt 



(t)-J2RnmklPkl(t) (10) 



kl 



where uj nm = (A„ -\ m )/K and max„ mfei \Re{R nmH )\ < 
min n ^ m |w nm |.The golden rule rates yield the following 
expression for the Redfield tensor 



Rn 



— 'Jim 1 nrrk ' u rik ^ 1 i rrrn 1 



-.(+) 
nmkl 



ikl ' 



(ii) 

With Ohmic dissipation with a Drude cutoff (w c ), the 
spectral function is 



Ji 



OLi TlLO 



1 + L0 2 /u) 2 ' 



(12) 



The formalism is valid provided a, < 1. Notice we have 
assumed that each qubit is individually coupled to an 
Ohmic bath. The rates are given by 



1 



coth( — — ) - 1 



i 



A 2 M ( ~ ) (w„ fe ,2) + A 1 Af 1)1 (13) 



where A* = o^\ m <J l z nk , being the matrix element 



of dg , and 



poo 

M ± (n,i) =V duj 
Jo 



uj 2 -n 2 



, , BTilo „ 
coth(- — ) ± 1 



(14) 



with V denoting the principal value. Similarly, 



-.(-) 



1 



Imnk 

i 



[A 1 Ji + A 2 J 2 



coth( 



A 2 M( + )( W;m ,2)+A 1 Af(+)( Wim ,l)l . (15) 
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Notice that the real parts of R contribute to the deco- 
herence and dissipation. For the purpose of studying the 
evolution of entanglement we will ignore the renormaliza- 
tion of the energy levels due to the imaginary part of the 
rates. This approximation has been verified to hold for a 
wide parameter regime for the Ising case 2 ^, and for weak 
coupling for single qubits coupled to Ohmic baths^i 

V. ENTANGLEMENT DYNAMICS 
A. J=0 



of this speed up can be traced back to the fact that the 
golden rule rates of the individual qubits are summed 
in the evolution equation for the density matrix and the 
factor of log(2) follows from the definition of entangle- 
ment. The fidelity of the state decays exponentially at 
a rate 21^ and the final state is a incoherent mixture of 
|t J,) and |4. T). Similar results are obtained if the initial 
state was the equivalent of the singlet or the Bell state 
ITT) + III)- Since we have dropped terms of higher order 
in n/e, this is not the complete story. For sufficiently 
long times the system reduces to a product state due to 
the presence of the bias e. 



Consider a system of two qubits that is prepared in 
an entangled state, e.g., the state with eigenvalue |J, 
which in the absence of A is the S z = triplet state. 
We will throughout this section use the terms "singlet" 
and "triplet" to refer to the states that are adiabatically 
connected to singlet and (S z = 0) triplet in the absence 
of A. Once initialized we allow the two qubits to evolve 
independently in the presence of dissipation, i.e., with 
J = and nonzero a. 

To leading order in the tunneling matrix element, we 
can write down the equations of motion for the density 
matrix. Given the initial condition, the only non zero 
elements are governed by the following equations, 



le an 

Pl3 = (y - al )pi 3 + — (P22 



P22 



-aTp 2 2 + aTp. 



33 



P33) 
: (P34 + P43) 



,ie . an. 
-(— + aT)p 13 + —{P22 



P31 = 

p'33 = aTp 2 2 - aTp 33 
The density matrix in the standard basis is 



P33J 






Vo 





_L 

2 

l p -2aK B Tt/h 
2 e 







-2aK B Tt/h 

_L 

2 




(16) 



(17) 



where terms up to order n/e have been retained. One 
immediately realizes that the time evolution of the den- 
sity matrix is controlled by the decoherence rate of the 
individual qubits, — a 5 , The dissipation is not 
present in the diagonal elements because they are not 
present at this order of approximation. The concurrence 
in this case is C — e -21 ^* and the entanglement is given 
by, 



E 



1 



e (-4iV) 1 
log 2 (- 



y/l -e(" 4I V) 



1 + V 7 ! - e(- 4r V) , 1 + v 7 ! - e(" 4r V) 

o lo §2 ( o 



)(18) 



The entanglement decays at a rate that is 2/ log(2) times 
faster than the decoherence rate of individual qubits as 
can be obtained by differentiating Eqn. 18. The origin 



B. Finite Coupling 

While the definition of entanglement can also be used 
for the case where the qubits remained coupled, the in- 
terpretation is less clear, as the two individual degrees of 
freedom are never truly isolated from each other. Here 
we rely on the observation that the exponential speed up 
in quantum algorithms is achieved due to the inability 
to reduce an entangled state to a product state. This 
notion can be applied even to the case where the qubits 
remain coupled. Here too we proceed by initializing in 
the same eigenstate as in the previous section. As seen 
in Fig. 1, for decreasing coupling the entanglement de- 
cays at a slower rate, but for larger coupling strength, the 
entanglement is completely lost and an incoherent state 
emerges at finite time. As the state evolves further, the 
entanglement becomes nonzero again. 

To get a better understanding of the dynamics, we con- 
sider the fidelity of the initial state (Fig. 2) . The fidelity 
characterizes the distance between two quantum states. 
It is defined as F =tr^p(Q) 1 / 2 p(t)p(Q) 1 / 2 . The entangled 
state evolves in a manner such that it becomes a mixture 
of the singlet and triplet states (note that the definition 
of singlet and triplet is only exact for A = 0). Ulti- 
mately the dissipation leads to the formation of product 
states with zero entanglement but, for low enough tem- 
peratures, the admixture with the singlet state is never- 
theless observed. For small coupling the state does not 
change greatly, as is reflected in the very slow decay of 
fidelity and entanglement. For larger values of J there is 
a drop in fidelity but it is rather featureless, as one might 
expect. Entanglement on the other hand displays a far 
more complicated behavior. 

The behavior of entanglement can easily be understood 
if one looks at the effect of dissipation on the quantum 
state. For the given Hamiltonian, the lowest energy state 
of the two qubits is indeed a product state. For short 
times, compared to the dissipation of the diagonal ele- 
ments of the density matrix elements, the "triplet state" 
has a lower energy compared to the "singlet" state. Thus 
the environment drives the system towards the singlet 
state leading to relaxation (i.e. loss of off diagonal ele- 
ments in the density matrix). The fidelity of the state 
tends to 1/2, which implies that the final state will be 
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FIG. 1: Entanglement as a function of time starting from 
the triplet-like state. Time is measured in units of t = h/e. 
The parameters used here are e = l,n = 0.2, T = 0.1 and 
a = 0.0001. 




FIG. 2: Fidelity as a function of time for the same parameters 
as in Fig. 



found to be in the initial state in half the measurements 
performed. Since the state evolves from one maximally 
entangled state to another, via an intermediate incoher- 
ent state, the entanglement reduces to zero before start- 
ing to increase. Note that for sufficiently long times, 
O(eto/n), the system will indeed lose all its entanglement 
content. 

To further illustrate the dynamics for long times, and 



FIG. 3: Entanglement as a function of time. Time is measured 
in units of t = h/e. The parameters used here are e = 1, n = 
0.2, J = 2 and a = 0.0001. 

the enhancement of dissipation with temperature, we 
now look at the time evolution of entanglement for J = 2. 
In Fig. 3, for low temperatures, the previous dependence 
is reproduced but for higher temperatures we notice that 
the long time decay increases rapidly while the initial loss 
of entanglement is hardly affected. Thus the intermedi- 
ate zero entanglement state is a reflection of the mixing 
of singlet and triplet state, whose rate is controlled by the 
coupling J, while the long time relaxation into product 
state is controlled by the temperature of operation. 

In all cases studied here, we have considered an initial 
state which is is adiabatically connected to the S z = 
triplet state for A = 0. One can study the effect of dissi- 
pation on the "singlet state" as well. Here the dynamics 
is rather trivial in that one only notices a decay to a 
product state. This is due to the fact that at finite J, 
there is no mixing with the triplet state as was observed 
before, since it has a higher energy. To verify that this 
physical picture is correct, we have checked that the roles 
of the states get switched when the sign of J is changed. 

C. Ising Coupling 

While interesting phenomena involving nonmonotonic 
behavior of entanglement is seen for the case where the 
coupling is isotropic, for an Ising-like coupling no such 
behavior is obtained (Fig. 4). We comment briefly on 
this case for comparison with the thorough study of gate 
fidelity under various operations in**. One might have 
speculated that the dissipation would lead to a quan- 
tum state with a lower energy, and hence even for the 
Ising coupling the nonmonotonic behavior persists, but 
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VI. CONCLUSIONS 




Time (t/t„) 

FIG. 4: Entanglement as a function of time for Ising coupling, 
starting from the entangled triplet state. Time is measured 
in units of t — h/e. The parameters used here are e = 1, n = 
0.2, K = 1 and a = 0.0001. 

this turns out not to be the case. The decay of entangle- 
ment to zero in Fig. 1 reflects the formation of the state 
(It J.) (til + lit) (lt|)/2. The density matrix for this state 
commutes with the Ising coupling term in the Hamilto- 
nian but does not do so for the SU{2) symmetric cou- 
pling. Thus, in the latter case, the state continues to 
evolve into an entangled state. While this argument is 
strictly valid at zero bias, e = 0, the computation per- 
formed here shows that it holds approximately even for 
finite e. The dynamics shown in figure 4 is representa- 
tive of all cases we have studied for difference values of 
the coupling K, both positive and negative. This is also 
consistent with the observation that the singlet is a pro- 
tected subspace for the Ising system, in that no mixing 
with the triplet is observedS*^ 



In this paper we have studied the time evolution of 
entanglement for two-qubit systems in the presence of a 
dissipative environment. In the absence of any coupling 
between the two qubits, entanglement of the initial state 
is found to decay in one limit at a rate that is about three 
times as fast as the decoherence rates of individual qubits, 
and numerical results suggest that this behavior is more 
general. In the presence of coupling the entanglement de- 
cay rate changes dramatically and depends on the nature 
of the coupling itself. For the case of an Ising like cou- 
pling between the two qubits, the loss of fidelity is faith- 
fully reproduced as a loss in entanglement. The same is 
not found to be true when the coupling between the two 
qubits is isotropic. For small coupling strengths, both 
entanglement and fidelity are lost rather slowly, but on 
increasing the the coupling not only does entanglement 
decay at a much faster rate, it also shows nonmonotonic 
behavior. In particular the entanglement is found to go to 
zero at intermediate times, before recovering on further 
time evolution. In agreement with the intuition that at 
higher temperatures and longer times, the system should 
evolve to an incoherent state, entanglement is found to 
vanish in this limit. 

This work is a first step in understanding the physics 
of coupled qubit systems in terms of their entanglement 
content. While the notion of entanglement of two qubits 
has been extensively studied and an analytic measure ex- 
ists, no such understanding for multiqubit systems exists. 
Such an understanding is vital for implementation of ef- 
ficient quantum algorithms. It is an important task for 
future investigation not only to quantify entanglement 
for more complex systems, but also to understand its 
evolution in the presence of a dissipative environment. 

The authors acknowledge helpful conversations with 
B. Plourde and T. Robertson, and support from DOE 
LDRD-366434 (V. A.), NSF DMR-0238760, and the Hell- 
man Faculty Fund (J. E. M.). 
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